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ABSTRACT 

L’vov and Falkovich (Physica D 57) have shown that the dimensionally possible inertial 
range scaling laws for Boussinesq turbulence, Kolmogorov and Bolgiano scaling, describe 
steady states with, respectively, constant flux of kinetic energy and of entropy. Following 
Woodruff (Phys. Fluids 6), these scaling laws are treated as similarity solutions of the 
direct interaction approximation for Boussinesq turbulence. The Kolmogorov scaling solu- 
tion corresponds to a weak perturbation by gravity of a state in which the temperature is 
a passive scalar but in which a source of temperature fluctuations exists. Using standard 
inertial range balances, the effective viscosity and conductivity, turbulent Prandtl number, 
and spectral scaling law constants are computed for Bolgiano scaling. 
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I. Introduction 


The spectra E u (k), E h (k), E 9 (k.) of buoyant turbulence axe defined so that 

roc 

MUi >= E u (k)dk 

Jo 

/•OO 

e 2 >= E e (k)dk 

Jo 


1 

2 <U ‘ 


< U36 > 


r °° 

/ E h (k)dk 
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where gravity acts in the 3 direction. The variables u, and 8 denote fluctuations about 
any mean velocity and temperature field. Dimensional analysis suggests two types of 
inertial range seeding for these spectra. If the gravitational coupling will be neglected, the 
temperature is a passive scalar and the velocity has the Kolmogorov spectrum 


E ~ e 2 / 3 k~ 5 / 3 


( 1 ) 


where e is the kinetic energy dissipation rate. Kolmogorov scaling can be expected to 
apply to wavenumbers k at which nonlinearity dominates the thermal forcing, so that 
Re(k) » Ra(k) for appropriately defined scale dependent Reynolds and Rayleigh num- 
bers. Although the limiting temperature distribution is uniform and the spectrum E e 
consequently vanishes in the absence of a source of temperature fluctuations, buoyant 
flows generally have such a source; in this case, there is a nonzero constant temperature 
variance dissipation rate N and dimensional analysis gives 

E 9 ~ N£~ 1/3 k~ 5/3 (2) 


In Kolmogorov scaling, the inertial range is isotropic, consequently the heat transfer spec- 
trum E h vanishes. 

Bolgiano identified a second possibility 1 in which velocity and temperature fluctu- 
ations are determined by N and g, where g denotes the product of thermal expansion 
coefficient and acceleration of gravity. Then dimensional analysis leads to 

E u ~ 0 4/5 iV 2/5 fc _11/5 

E h ~ g 1/5 N 3/5 k~ 9/5 (3) 

E 9 ~ g~ 2 / 5 N*/ 5 k~ 7 / 5 
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These scaling laws remained largely a dimensional possibility until L’vov and Falkovich 2 
clarified their possible dynamic significance: whereas Kolmogorov scaling corresponds to a 
steady state inertial range with constant energy flux e, Bolgiano scaling corresponds to a 
steady state inertial range with constant entropy flux JV; the identification of temperature 
variance dissipation with entropy flux is due to L’vov . 3 This scaling applies when thermal 
forcing dominates the nonlinearity so that Ra(k) » Re(k). A nonzero heat transfer spec- 
trum is possible for Bolgiano scaling, in which the force of gravity introduces a preferred 
direction. This physical picture might suggest that Kolmogorov scaling will be observed 
in forced convection and Bolgiano scaling will be observed in free convection; however, 
arguments 3 that measurements in very high Rayleigh number Rayleigh- Benard convection 
experiments are consistent with Bolgiano scaling cannot be considered conclusive . 4 Nev- 
ertheless, turbulence sustained even at small scales by buoyant forcing is an interesting 
theoretical possibility which deserves investigation. 

This paper considers Kolmogorov and Bolgiano scaling following Woodruff 5 as simi- 
larity solutions of the direct interaction approximation (DIA) for buoyant turbulence . 6 It is 
shown that Kolmogorov scaling can be treated as a perturbation of the passive scalar state 
in which the gravitational coupling g vanishes. For Bolgiano scaling, this viewpoint results 
in close connections with the e-expansion of Yakhot and Orszag , 7 which Woodruff inter- 
prets as arising from an asymptotic evaluation of the integrated DIA response equation. In 
Eulerian theories, the e-expansion also serves as a scheme of infrared regularization, which 
is required since any attempt to compute amplitudes in the spectral scaling laws using 
the DIA equations directly is defeated by the well-known infrared divergence . 8 The only 
fully satisfactory solution of this divergence is the formulation of a Lagrangian theory . 9,10 
However, the construction of analytical solutions for the Eulerian theory requires some 
substantial approximations; the added complexity of a Lagrangian theory of any problem 
with coupled fluctuating fields may justify more or less ad hoc modifications of the Eulerian 
theory like the e-expansion. 

II. Direct interaction approximation for Boussinesq turbulence 

It will be convenient to write the Boussinesq equations in matrix form as 

G 0 -‘U(t) = r(i) f dpdqV(p)V(q) (4) 

J k=p+q 
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In Eq. (4), U is the vector 



Written explicitly, the nonlinear term on the right side of Eq. (4) is 


where 6 


In Eqs. (4)-(7), 


r imn (k) dp dq U m (p)U n (q ) 

J k=p+q 


1 

2 


Pimn{ k) if i, m,n^4 
k m if i = n — 4 or i = m = 4 

0 otherwise 


A = (ft,k) p = («,p) 5 = (J2-w,q) 


(5) 

(6) 
(7) 


and 

Pimni k) = fc m P in (k) + k n P im ( k) 
P,y(k) = 8ij - kikjk~ 2 
The matrix G 0 in Eq. (4) is defined by 


Go(k) 1 = 


(-iQ + u 0 k 2 )l gP 3 

0 (—iCl + Kok 2 ) 


(8) 


where the vector P3 has the components P l3 and uq and kq are the molecular viscosity and 
conductivity respectively. Inverting the matrix in Eq. (8) gives the bare Green's matrix 


Go {k) — 

where Gq and Gq are the bare propagators 


GS(k) I -G5(k)G*(k) ff P 3 

0 Gg(i) 


(9) 


GS(i) = (-in + Uok 2 )- 1 

G e 0 (k) = (-in + Kok 2 )- 1 (10) 


It should be noted that, unlike the passive scalar equation in which g = 0, the Boussinesq 
equations are effectively nonlinear in T. 

The direct interaction approximation (DIA) for Boussinesq turbulence has been de- 
rived by Kraichnan. 6 We will attempt to construct an approximate analytical solution to 
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the DIA equations in the inertial range. The descriptors of buoyant turbulence in the DIA 
are the correlation functions Q±j, Qf , Q T defined by 

< ui{k)u } {k') > = Q u tJ (k)8{k + k') 

<Ui(k)T(k')> =Q}(k)S(k + k') 

<T(k)T(k')> =Q T (k)S(k + k') (11) 

and the response or Green’s functions, 

<?&(*) =< Sui(k)/6ff(k) > 

Gf(k ) =< Sui(k) / 6 f T (k) > 

Gj u (k) =< ST(k)/Sff(k ) > 

G e (k) =< 8T{k)l8f T (k) > 

where 8ff and 8f T denote small perturbations added to the velocity and temperature 
equations respectively. 

The analysis will be based on the Langevin equation representation of DIA which for 
an inviscid steady state takes the form 

-iSlui(k) + gT{k)P a { k) = -r £(*)u r (*) - rjfCk)T(k) + /“(*) 

-iQT(k) = —r}J u (k)u r (k) - r] T {k)T{k) + f T {k) ( 12 ) 

The damping factors 77 are defined by 

rll(k) = jiWk) / <ip dq r dm {P„„(p)G*(p)0;.(«) + P.,.(q)G*(fl%.(#) 

^ J k=p+q J — 00 

- PrG'J (P)Q h M) ~ qrGl T (p)Ql(p)} ( 13 ) 

n‘ T (h = - jPimn(k) [ dp dq I dm \p,G u J{p)Q‘Jq) 

^ v/k=p+q J — oc 

+ «rG; r (s)« m (p)} (14) 

’lj x ‘(k) ~ \ j dp dq I dm \k„p,G , (p)Q''Jq) + k„q r G* (q)Q*(p) 

+ k m P mr ,(p)G‘(p)Q*(q) + k„P mT ,(q)G*(q)Q k ,(p) 

- k„Gl T (p)Q T (q) - k n G° T (q)Q T (p) ( 15 ) 

- k m P„ r ,(q)Gl“(p)Q' m ,(p) - k„P nr ,(p)Gl“(q)Q* m ,(q) 

- k m q r G^(p)Q‘ m ,(p) ~ kmPrGl-miM)} 
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( 16 ) 


rj T (k) = j f dp dq f du> {k n p m G 6 (p)(?ro n (g) + k n q m G 9 (q)Qmn(P) 

4 J k=p+q J — oc 

- k nPr G u n T (p)Q h r (q ) - k n q r G u J(q)Q h r (p)} 

In Eq. (12), /“ and f T are random forces with the correlations 

< /*«/“(-*) > = F““(i) 

< rn'm-h > = f“ T (i) 

< T(k)T(-k) > = f TT (t) 

where 

F TT (k) = k m k r j dp dq Q u mr (p)Q T \q) + Qt(P)Qr(q) 

J k=p+q 

F t uT (k ) = k n P irs ( k) / dp dq Q h r (p)Q u ns (q ) (17) 

= P imn (k)Pj rs (k) / dp dg Q^ r (p)Q“ a (9) 

Jk=p-\-q 

The DIA theory is completed by the definitions Eq. (11) of the correlations in terms 
of the fields u , T and by the equations of motion for the G’s, which are given by 

-iS2G? ; (i) + gGj‘(k)P a (k) + ,« r Gj“(t) + (i)G^ = 1 

-iiiGf(k) + 9 G r (ifc)P i3 (k) + gf(k)G T (k) + tf r (k)Gf(k) = 0 
-iSiGp(k) + r, T (k)G]'(k) + r,J‘(k)G‘ TJ Ck) = 0 
-iflG'fi) + ^(iciG’Ck) + gp(k)Gf(k) = 1 ( 18 ) 


These equations are simply the conditions that the solution of the Langevin equations, Eq. 
(12) is given by 


f u ( k) 

U T (k)\ 


(19) 


U(Jfc) = G (k) 

in the matrix notation of Eqs. (4)-{9). In Eq. (19), G is the renormalized response matrix 

G(fc) = 

The notation means that G r “ is a row and G uT is a column. 


G,,(k) G° T Ck) 
lGp(i) G«(i) 1 
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The DIA equations are considerably more complicated than the baxe equations; al- 
though the coupling r?“ r can be considered a renormalized gravitational coupling, rj Tu has 
no analog in the bare equation. Moreover, whereas the bare theory contains only the two 
response functions of Eq. (10), DIA adds the functions G uT and G Tu . Analytical progress 
will clearly require rational simplification of these equations. It will be convenient to solve 
the DIA equations by an iterative procedure, starting with an ansatz for the G's and F’s. 
When these are known, the correlations can be computed from Eqs. (11) and (19). Know- 
ing G and Q, substitute in Eqs. (13)-(16) to obtain the damping factors. This procedure 
can now be iterated, since an updated G can be found by solving Eq. (18), and updated 
forces can be computed from Eqs. (17). 

The initial assumption for the Green’s matrix will be 


G(k) = 


G u (k) I -G u (k)G e Ck)gP 3 
0 G 9 (k ) 


( 20 ) 


so that the bare and renormalized Green’s matrices have the same structure, and are 
determined by two scalar Green’s functions G u (k) and G 9 (k) alone. This type of diagonal 
approximation is useful in treating coupled field problems. Leaving the random force 
general temporarily, the correlation matrix will be determined by forming correlations of 
the amplitudes found from Eqs. (11) and (19). The explicit expressions are 


Q0) = G*Ck)G u (-k)F t y(k) 

- gG u Ck)G u (-k)G e (-k)PM^)F t uT (k) 

- gG u (-k)G u (k)G 9 (k)P iz (k)Ff(k) 

+ g 2 G u (k)G e (k)G u (-k)G e (-k)Pi 3 (k)Pj 3 (k)F TT {k) 

Qi{k) = G u (k)G 9 (-k)F* T (k) - gG u Ck)G 9 (k)G 9 (-k)P l3 (k)F TT Ck) 

Q T (k) = G e (k)G 9 (—k)F TT (k) (21) 


Now we will adopt for Kolmogorov scaling the ansatz 

Fjj u (k) = Pij(k)Cpek~ 3 F uT — F TT = 0 (22) 

and for Bolgiano scaling 

F TT {k ) = CpNk~ 3 F uu = F uT = 0 (23) 
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The forces in Eqs. (22) and (23) are white noise in time. It will be shown later that a 
corrected force can be derived which is not white noise in time. Thus, DIA predicts that 
the effective forcing is in fact certainly not white noise in time; the introduction of such 
a forcing is merely used to provide an initial approximation to the solution of the DIA 
equations in the inertial range. 

A. Kolmogorov scaling 

Substitution of Eq. (22) in Eq. (21) results in 

Qij(k) = G u (k)G u (—k)Pij(k)Cpek~ 3 

Qi(k) = 0 (24) 

Q T (k ) = 0 

These are the correlations appropriate to a purely passive scalar. The vanishing of the 
Q T correlation reflects the absence of a source for temperature fluctuations: the long time 
temperature distribution is uniform and therefore vanishingly small. 

Substituting Eq. (24) in Eq. (13) leads to the isotropic renormalized damping 

»$(*) = 

where 

!,“(*) = jP.mn(k) / ft.r.(p)*..(q) /” G*(p)G*( $)G* (-«)C? Bq' 3 (25) 

^ J k=p+q J — oc 

Substituting Eq. (24) in Eq. (14) leads to the renormalized diffusivity 

P.m(q) r G'(p)G*( 9 -)G-(-i)Cfet-* (26) 

^ »/k=p-hq J —oc 

Substituting Eq. (24) in Eq. (11) leads to the renormalized gravitational coupling 

,“ r (i) = ,* r (t)P a (k) 

where 

V UT (k) = -\gP 3 mn(k) I Pr-Pm3(p)P r „(q)x 

4 J k=p+q 

/ oo 

G*(p)G*(p)G“(9)G“(-i)Cj£ 9 - 3 (27) 

■ OO 
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The coupling r) Tu vanishes identically. 

In this approximation which neglects G Tu and G uT , the system of equations Eq. (18) 
for the response functions is replaced by the simplified system 

-ittG u (k) + rj u (k)G u (k) = 1 

-inG e (k) + r 1 T Ck)G e {k) = l (28) 

The couplings rj uT and rj Tu do not appear in these simplified response equations, and can 
therefore be ignored. 

It is shown in the Appendix that the assumption of white noise in time forcing is 
consistent with Markovian damping in which 

7 ? “(fc)=7 ? “(fc) 
rj T (k) = r/ T (k) 

The corresponding response functions are exponential: 

G u (k) = (-in + ^k))- 1 

G 9 (k) = (—iQ, + ar/(A:)) _1 (29) 

The damping factors ri(k ) and r] T (k) for Kolmogorov scaling are 

77(fc) = C D e 1/3 k 2/3 

i f(k) = aC D e 1/3 k 2 ' 3 (30) 

The constant C D is to be calculated from the theory. Eq. (30) assumes the existence of a 
scale independent turbulent Prandtl number 1/a, which is also to be computed from the 
theory. Substitution of Eqs. (29) and (30) in Eqs. (25) and (16), evaluating the integrals 
by the e-expansion following Woodruff 6 leads to the results for the passive scalar previously 
obtained by Yakhot and Orszag. 7 In view of the vanishing of Q T and Q h , Eq. (17) shows 
that the ass um ption Eq. (22) for the random force was consistent with DIA. 

The passive scalar solution obtained in Eqs. (25)-(26) can be corrected by adding a 
thermal force to Eq. (22): 

F% u (k) = Pij(k)C£ek~ 3 F uT = 0 F TT = C^Nk~ 3 
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and treating g as an expansion parameter. Substituting this forcing in Eq. (21) and 
integrating over u> leads to 0(g 2 ) corrections to the correlation Q u and to the results 


Q k 

q t 






N 

? £ 2/3 

JL* 

- 1/3 K 


k-y* 

- 5/3 


The result for Q T is equivalent to Eq. (2). In this analysis, nonvanishing coupling to the 
velocity equation is needed to produce a nonvanishing temperature spectrum. 


B. Bolgiano scaling 

Substituting Eq. (17) in Eq. (15) leads for Bolgiano scaling to the correlations 


Q-.(k) = g 2 G u (k)G\-k)G 9 Ck)G 9 (-k)P t3 (k)P j3 (k)C^Nk- z 
Q h t {k) = -9G u (k)G 9 (k)G 6 (-k)P i3 (k)CjNk- 3 

Q T (k ) = G 9 (k)G 9 (-k)C^Nk~ 3 (31) 


Eqs. (31) imply various relations among inertial range spectra which can be derived by 
integrating over frequency to obtain single time correlations. Again assuming Markovian 
forms Eq. (29) for the response functions G and G e , the single time correlations are found 
to be 


/ °° - C T Nk ~ 3 

dw Q T (k ) = Cf 


-OO 

roo 


2arj(k) 

/ OO £» 

duQ h (k) = -gP l3 (k) — 

P i3 (k)Q T (k) 


C^Nk-* 

+ l)v(k) 2 


(a + l)rj(k) 

/ °° 

dw Q T (k ) = 

“OO 


(a + l)rf(k) 2 


2 a(a + 1) rj(k) 3 
P, 3 (k)P j3 (k)Q T (k) 


CpNk~ z „ „ „ , 

F P l3 (k)P ; - 3 (k) 


(32) 


Thus, all of the single time correlations can be expressed in terms of Q r (k). Substituting 
these correlation and response functions in Eqs. (13)-(16), we obtain the renormalized 
quantities 


ri(k) = \g 2 Pimn{k) f {prPm3{p)Pn3((l)CpNq 3 X 

4 J k=p+q 
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/ oo 

du G*(p)G 6 (p)G*(q)G e (q)G*(-q) 

-OO 

/ OO 

<^G“(p)G“(<,)G“(-<j)G'(<i)G 4 (- ? -)}( 33) 

-oo 


r li T {ty — ^ gPimn 


(k) f p r Pm 3 (p)Pr 3 (q)Pnz(q)CfNq- :i X 

J k=p+q 


/ oo 

d w G“(p)G , (p)G’(«G“(-«)G 4 (5)G 4 (-?) 

-oo 


(34) 


q T (k) = j g 2 k « p m P m i l q)P 7i i<q)CfNq 3 / dui G*(p)G u (q)G u (-q)G*(q)G e (-q) 

J k=p+q J — oo 

/ oo 

<t> G“(p)G 4 (p)G“( 9 -)G 4 (,j)G 4 (-<;) (35) 

-oo 


.£*(*) = I 9 K Pm P„ 2 (q)C T F N q - 3 r du G\p)G'‘(q)G\q)G‘(-q) 

J k=p+q J — oo 

/ oo 

<tp G‘(p)G'(q)G s (q)G e (-q) 

-oo 

/ CO 

du G'(p)G%p)G‘(p)G l (q)G > (-q) (36) 

-oo 

As in the case of Kolmogorov scaling, the couplings fj 1 ” and p uT can be ignored at this 
level of approximation. 

III. Determination of inertial range constants for Bolgiano scaling 

A. The turbulent Prandtl number 

The procedure for determining inertial range constants is standard. 11 Because of Eq. 
(32), there are only three inertial range constants, namely the constant Cd in the time 
scale defined for Bolgiano scaling by 


Tf(k) = C D g 2/5 N^ 5 k 2 ^ 


(37) 


the turbulent Prandtl number 1/a defined so that the thermal damping function r/ T (k ) is 


t 7 r (fc) = aC D g 2/5 N 1/5 k 2 ' 5 


(38) 
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and the constant C'g in the spectral scaling law 


E r {k) = C B g~ 2/5 N 4/5 k- 7/5 (39) 

Setting ft = 0 in the approximate response equations Eq. (28), 

G K (k) v H (k) |n=o= 1 

G 9 (k)rj T (k ) |n=o= 1 (40) 


Performing the frequency integrations in Eq. (40), 


1 = 


4 7](k) 


»»( k ) / 

J k=i 


dpdqx 

= p+q 

arj(p) + (2 a 2 + 2a + 1 )rj(q) 


PrPm3(p)Pn3(<l)Q T {q) 


a(a + 1 )n(q)(v{p) + q{q))(v{q) + ar/(p))(r)(p) + ar](q)) 


+ 


(a + 1 Hq) + r]{p) 


(a + 1 )v{q) 2 (v(p) + p(q))( ar i(q) + vip)) 


Pmrs(p)Pn3(q)Ps3(q)Q T (q) 


(41) 


1 = 


L. 


dpdqx 


knPmPm3(ci)Pn3(q}Q T {q) 


4cn?( k ) jk=p+q 

<Mp) + (a + 1 )v(q) 

a(a + 1 )v{q) 2 (ar](p) + Tj(q))(i ](p) + rj{q )) 

+ [a(a + i)v{q)(p(p) + q(q))(an(p) + p (q))(v(p) + Mq))) 

knPrPnz(p)Prz{q)Q T (q) 


-1 


X 


(42) 


Substitution in Eqs. (41) and (42) of the inertial range forms Eqs. (3) leads to the 
usual problem that the integrals are infrared divergent. We follow Woodruff’s treatment of 
this problem for hydrodynamic turbulence. 5 We begin by seeking scaling forms for E e and 
■q such that the right hand sides of Eqs. (41)-(42) are ultraviolet divergent. The scalings 
are found by assuming that 

rj ~ & 2-f / 5 

The right sides of Eqs. (41) and (42) must be scale independent. Thus, if Q(k ) ~ k A , 
necessarily A = 3 — 4e/5 so that 

r)(k) ~ fc 2-f / 5 

E°(k) ~ fc 5_4f/5 (43) 
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Bolgiano scaling, defined by Eqs. (3), (37), and (38) corresponds to e = 8. Substituting 
the scaling forms of Eq. (43) in Eqs. (41) and (42) shows that the integrals are ultravi- 
olet divergent for e < 0 and logarithmic when e = 0. The values of the integrals when 
e > 0 can be approximated by asymptotic expansion about e = 0. Because of the ultravi- 
olet divergence, the integrals are dominated in this case by distant interactions in which 
k/p,k/q — *• 0 and p, q — * oo. Substituting the scaling forms Eqs. (43) in Eqs. (41) and 
(42), evaluating the integrals in the distant interaction limit, and performing the spherical 
integrations leads after lengthy but straightforward calculations to the results 

ff 2 p ( k \ f* ot ?(p) + ( 2q2 + 2ct + 1 ) ? ?(g) 

4 r}{k) rmn J a(a + l)q{q)(q(p) + r}(q))(rj(q) + ar}(p))(r)(p) + ar](q)) 
PrPni3{p)Pn!l( , i)Q (?) 

= ^§f [ “ 2( “ + 1)Sl “ 1 i (6 ° 2+5 “ + 3) l (44) 

g 2 p (k] f* (« + i)g(g) + g(p) 

4 q(k) rmn J (a + l)rj(q) 2 (r](p) + rj(q))(ar](q) + ??(p)) 

Pmrs(p)Pn3(q)Ps3(<i)Q T {q) 

= 16 + 15^ 6a2 + 18a + 16 ^6e ( 45 ) 


g 2 p (k) [* «g(p) + (q + iMg) 

4 ar](k) rmn ’J a{a + l)q(q) 2 (ar](p) + p(q))(v(p) + p(q)) 

knPmP m3 (Q)f”n3 (q)Q (g) 

= I5^ l<lJ(0 + 1)2| ' , I (2a + 1) l (46) 


g 2 p (k] r Q7 ?(p) + (2a 2 + 2a -f l)p(g) 

4 r}{k) rmn J a{a + l)q(q)(rj(p) + r/(q))(ari(p) + T](q))(r}(p) + aq(q)) 
knPrPn3{p)Pr3(<i)Q T (q) 


( 47 ) 


where 


r=| t / *>+/ 

J £ Jq=k — p,p>fc J p 


dq] 
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Note that the symmetrization in f* is required because the original integration over 
wavevector triads is certainly symmetric in p and q, but J^_ k > fc dp by itself for exam- 
ple is not symmetric. The restrictions p > k and q > k in f* are infrared regularizations 
which follow from the e-expansion. 5 This regularization suppresses a crossover to a domi- 
nant sweeping interaction when e = 7 which reflects the infrared divergence characteristic 
of Eulerian theories. 6 Evaluation of the integrals can be simplified by summing over the 
direction of gravity, which always enters these integrals in pairs. 

Substituting the results Eqs. (44)-(47) in the integrated response equations Eqs. (41) 
and (42) leads finally to 


= i(a 2 (a + l) 3 ]-‘|- ±(12a 3 + 48a 2 + 42a + 6) 
C b 16 be 15 


£t = l[ a 3( a + l)]-il?4 

C B 16 6e 3 


(48) 


Equating the right hand sides of Eqs. (48) leads to an equation in a alone. Clearing 
fractions in this equation, 


12a 4 + 48a 3 + 2 a 2 - 74a - 40 = 0 


which has a unique positive root 

a ~ 1.275 

Now substituting in Eq. (48), we find 

^ = A = 0.0037 
Cb 


(49) 


B. Evaluation of the constants Co an d Cb 

Following Kraichnan, 12 a second relation between Cq and Cb is found from the con- 
stancy of inertial range entropy flux defined 3 by 

N = f dk f dp dq S(k, p, q) — / dk f dp dq S(k, p, q) (50) 

Jk<ko o J 2ko>k>k 0 

where the entropy flux density 5 is 


S(k,p,q) = -ik m < ti m (p,f)T(q,f)T(-k,<) > 
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The flux N should be finite and independent of k 0 in the inertial range. In the DIA closure, 


S(k,p,q) = lim<_ooS(k, p, q, t) 


where, reverting to time domain notation for simplicity, 

E(k,p,q,<) = ~^k m j ds {G“(p,#,s)P mrs (p)Q*(q,s,t)Q^(k,s,t) 

+ G 0 (ci,t,s)q r [g^(p,t,s)<5*(k,5,t) + g^ r (p,t,s)g T (k,t,s)] 

- c <> (k,t,s)fc r [g^(p,f,s)g|:(q,5,<)] + g" r (p,*,s)g T (q,M)] } 


+ 


ds jf 


drx 


{ PrG d (p , 5 , r ) [Q* (q, r, t)Q T { k, r, t ) + g fe (k, r, t)Q T (q, r, t)] } 

+ \gk m I ds j drx 
* Jo Jo 

+ {<lrG 9 (q, t , r) [Q*(p, r, s)g T (k, r, t ) + (k, r , t)Q T ( p, r, s)] 

- k r G e {k , f, 5 ) [Q* (p, r, s)g T (q, r, t) + Qr (q, r, f )Q T (p, r, s)] } 

The time integrations are performed using fluctuation dissipation relations, which follow 
from Eqs. (29) and (19) 

< ,,(k, (1 )T(k >= «(k + ko { £ *4 > ;; 


Then a lengthy but straightforward calculation leads to 


*5(k,p,q) = -^©l(fc,P,?){ fcmi 3 mr S (p)-Pr3(q)^3(k)g ,l (A;)g fc (9) 

+ Pmr(p)^ m g U (p) [9rQ T (^) “ KQ T (q)] 

+ Pml(p)kmQ h (p)[qrPrz(k)Q h (k) - k r P r3 (q)q h (q)] } 

+ 7 ;gkmPmz(p){®i(k,p,q)& 2 {k,p,q)[prPrz(q)Q h (q)Q T (k) + PrPrz('k)Q h (k)Q T (q)] 
+ [&i(k,p,q) 03 {k,p,q) + (2T](p))~ 1 0i{k,p,q)}q r P r3 (p)Q h (p)Q T (k) 

+ [0i(fc,p,«)02(*,p,«) + ((1 + a)^(p)) _1 0i(fc,P,9)]9rPr3(k)g fc (fc)g r (p) 

- [01 {k,p,q)Q 2 (k,p,q) + ( 2ri(p))~ 1 Q 1 (k,p,q)]k r P r3 {p)Q h (p)Q T (q ) 

- [0i(fc,p,g)0 2 (fc-P,?) + ((1 +a)ri(p))~ 1 ®i{k,p,q)]k r P r3 (q)Q h (q)Q T (p)} (51) 
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where the time scales O z (k,p, q) are defined by 


01 (k,p,q) = [r](p) + ar](q) + ari(k )}- 1 

©2 (k,p,q) = a~ 1 [r]{p) + T](q) + rjik )^ 1 

03 (k,p,q) = [ar](p) + r)(q) + aq(k) ] _1 (52) 


Note that the time scales in the g dependent term in Eq. (51) have a nonstandard form, 
which might not suggest itself in a phenomenological approach like EDQNM. 

Now Eqs. (37)-(39) can be substituted in Eq. (51), with Eq. (32) used to express Q h 
and Q u in terms of Q T alone. To evaluate the flux integral Eq. (50), define the quantity 
T(k,p,q ) by 

/ dp dq T(k,p,q) = j dp dq S(k,p,q) 

J J k=p-fq 

The scalings are such that 


T(\k,\p, \q) = \- 3 T(k : p,q ) 

a result equivalent to the choice of a k~ 3 Langevin force in Eq. (23). This force scaling is 
also equivalent to Orszag’s condition 13 on the scaling laws for 0 and Q that a stationary 
inertial range is possible. Accordingly, the -3 forcing appears to be universal, a suggestion 
made on other grounds in Ref. 14. This scaling also implies that the inertial range transfer 
does not depend on the wavenumber ko and permits the transfer integral to be reduced to 
a double integral over a finite region following Kraichnan. 12 Numerical integration of Eqs. 
(21)-(22) leads to 

= B = 0.21 (53) 

Eqs. (49) and (53) imply tentative values of the inertial range constants 

C B = A 3 / 5 B i/5 = 0.01 
C D = A 2 / 5 B l/5 = 0.08 

Whereas the turbulent Prandtl number 1/a = 0.78 is not far from the passive scalar value, 
these constants depart markedly from values for Kolmogorov scaling. Further investigation 
of the inertial range constants is indicated. 
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IV. Iterative solution of DIA 

The approximation considered so far replaced the exact system of equations for the 
response functions Eq. (18) by the simplified set Eq. (28) by ignoring the response 
functions G Tu and G uT . This approximation also entails ignoring the couplings functions 
r{ Tu and ry uT . But Eqs. (13)-(16) provide expressions for rj Tu and tj uT . As noted previously, 
these can be substituted in Eqs. (18) to compute corrected response functions including 
G Tu and G uT . The force correlations can be corrected similarly. It is evident that this 
iteration will not alter the scalings, but changes to the computed inertial range constants 
cannot be ruled out in advance. 

To indicate the general character of this iterative procedure more explicitly, it will be 
convenient to simplify the spatial dependence of the structure of Eq. (18). Assume that 
the solution has the form 

<?&(*) = Gi&)Pii(- k) + G 2 “(W 3 (k)P > 3 (k) 

Gf(k) = G' T (k)P e ,(k) 

Gj'(k) = G T "(k)P 3 j (k) (54) 

The results from Eqs. (33)-(36) have the corresponding forms 

= r,W)Pij( k) + -J 2 “(*)P«(k)P > 3 (k) 

V? T (k) = ,“ 3 '(k)P is (k) 

t ,f(i) = <l T "(k)P 3j ( k) (55) 

where the Markovian character of the approximation has been used. Now substitute Eqs. 
(54) and (55) in Eq. (18) and note, as observed previously, that the index 3 corresponding 
to the direction of gravity, always appears in obvious pairs. The isotropic part of the 
solution can be evaluated by summing over this index. 15 This procedure leads to a system 
of linear scalar equations with a simple structure 

[~iQ + rj u {k)]G u {k) + rj uT (k)G Tu (k) = 1 
[-in + rj u (k)}G uT {k) + r/ uT (fc)(? r (]fc) = 0 
[-in + T) T (k)]G Tu (k) + T} Tu (k)G u (k ) = 0 
[-in + 7 ? T (it)]G r (it) + V Tu (k)G uT (k) = 1 (56) 
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where G u = G“ + G% and 77 “ = 77 “ + 772 . The notation in Eq. (56) has been simplified by 
including the bare gravitational coupling in r] uT . Eq. (56) can be solved for the response 
functions: 


G u (k ) = [-?fl + 77 T ]/[(-7'fl + T] u )(-iQ + 7] T ) - 77 uT ? 7 Tu ] 

G uT (k) = -ri uT l[{-iV + + r} T ) ~ V“] 

G Tu (k ) = -rj Tu /[(-iQ + 77 «)(-*fi + rj T ) - 77 " T 77 r “] 

G*(k) = [-til + ri U }/[(~in + r} U )(~M + V T ) ~ V uT rj Tu ] (57) 

where the damping functions 77 are found by substituting the lowest order response func- 
tions Eq. (29) in Eqs. (33)-(36) and setting Q = 0. Whereas the lowest order response 
functions are rational functions of first order, the corrected response functions of Eq. (57) 
are rational functions of second order with poles at the points 

n = -i{ 77“ + T7 r ± [(77“ - 77 t ) 2 + 477“ V“] 1/2 } 

These response functions can be substituted in Eqs. (33)-(36) to produce a system of 
algebraic equations for corrected damping functions 77 . Continued iteration between Eq. 
(57) and Eqs. (33)-(36) produces a sequence of approximations to both the G’s and the 
77 ’s of increasing complexity. This procedure is loosely analogous to the representation of 
correlation functions by continued fractions , 16 but is evidently much less explicit. 

Eq. (57) also suggests an alternate initial approximation for the response matrix 
to Eq. (20). Namely, substitute Eqs. (37) and (38) for rj u and r] T , and introduce the 
remaining damping functions 

rj uT (k) = C u D T g 

r] Tu (l t) = Cl u g- 1/5 N 2/5 k~ 6 ^ (58) 

Evaluating Eqs. (18) when Q, = 0 will produce a system of equations for all of the inertial 
range constants analogous to the system derived from Eq. (40). The constants should be 
evaluated more accurately by this procedure, but closer investigation must be left to future 
research. From the viewpoint of the e-expansion, it is easily verified that rj uT is 0(e) and 
that rj Tu is irrelevant for e < 5. 
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In outlining the higher order approximations, we have adhered to the distant interac- 
tion limit with its attendant Markovian damping. An interesting alternative is to drop this 
limit while retaining some kind of infrared regularization. Thus, substituting exponential 
response functions in Eqs. (13)-(16) makes these damping functions non-Markovian. The 
response functions computed from Eq. (18) with these damping functions will no longer be 
exponential. Similarly, the corrected forcing is not white in time. Like the iterative solution 
outlined above, this possibility of correcting the distant interaction limit and evaluating 
the correct non-Markovian damping is an interesting subject for future investigations. 

V. Conclusions Bolgiano and Kolmogorov scaling have been treated as solutions of the 
direct interaction equations for the inertial range of buoyant turbulence. A preliminary 
computation of the universal inertial range constants for Bolgiano scaling has been com- 
pleted. The solution entails some simplifications of the time dependence of correlation 
functions and response functions in the inertial range, namely that the effective forcing is 
white noise in time and that the effective damping is Markovian. An iterative scheme which 
treats the present solution as a first approximation and which corrects these hypotheses is 
proposed. 

APPENDIX Markovianized DIA 

Let the real function If(£), 0 < £ < oo satisfy 

pOO 

H( 0) = 1, H( 0 < 1 for £ > 0, / < oo (54) 

Jo 

Then standard properties of delta functions imply 

XH(X(t-s)) ~ S(t - s) for A -► oo (55) 

Rewrite Eq. (25) in the time domain, and evaluate the wavevector integrals in the distant 
interaction approximation in which k — *• 0, p, q — ► oo. Then 

p u (M,s)= f dpdqB(k,p,q)G u (p,t,s)Q u (q,t,s) 

J k=p+q 

/ * ^p 

dp { km -Q— (k, p, p)G“(p, t, s)Q u (p, t, s) 

dQ u 

- B(k,p,p)G u {pJ,s)k m p m p 1 -^-(p,M)} 
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where B(k,p,q) denotes the product of projection operators in Eq. (25) Assuming time 
stationary similarity forms G u (p, t, s) = G u (p r {t — s )) and Q u {p, t, s) — R u (p r (t — s))Q u (p) , 
the properties Eq. (54) of H may reasonably be postulated of the product G U R U . There- 
fore, Eq. (54) implies that in this limit the damping is Markovian 

r/ u (M,s) = 6(t-s)i 7 u (k) 

and Eq. (28) implies that the Green’s function is exponential, 

G u (k,t,$) = exp [(s - f)r?“(k)] for t > s (55) 

Likewise evaluating the force correlation Eq. (17) in the distant interaction limit implies 
that the forcing is white noise in time: 

<f?(k,t)f?(k',s)> = 6(t — s)£(k + k')F“j(k) (56) 

Computing the correlation function from the relation 

t ft 

dr i G u (k,t , ri ) / dr 2 G u (k',t,r 2 )x 

Jo 

<r(k,n)f(W,r 2 )> 

using Eqs. (68), (69) shows that the fluctuation dissipation relation 

Q u (k,t,s) = Q u (k)[G u (k 1 t,s ) + G u (k,s : t)] 

expressing the time dependence of the correlation functions in terms of the response func- 
tion is also valid in this limit. 
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